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ABSTRACT 

In this Letter we investigate the properties of advection-dominated accretion flows (ADAFs) 
fed by the evaporation of a Shakura-Sunyaev accretion disk (SSD). In our picture the ADAF 
fills the central cavity evacuated by the SSD and extends beyond the transition radius into a 
coronal region. We find that, because of global angular momentum conservation, a significant 
fraction of the hot gas flows away from the black hole forming a transsonic wind, unless the 
injection rate depends only weakly on radius (if r 2 a oc r - ^, £ < 1/2). The Bernoulli number 
of the inflowing gas is negative if the transition radius is <^ 100 Schwarzschild radii, so matter 
falling into the hole is gravitationally bound. The ratio of inflowing to outflowing mass is ~ 1/2, 
so in these solutions the accretion rate is of the same order as in standard ADAFs and much 
larger than in advection-dominated inflow/outflow models (ADIOS). The possible relevance of 
evaporation-fed solutions to accretion flows in black hole X-ray binaries is briefly discussed. 

Subject headings: accretion, accretion disks — black hole physics — hydrodynamics 



1. INTRODUCTION 

In recent years X-ray and optical observations pro- 
vided increasing evidence that advection-dominated 



flows (ADAFs, Ichimaru |1977| ; Narayan k Yi |1995b| ; 
Abramowicz et al. 1995 ) are a ubiquitous feature of 
accretion at all scales, from black holes (and possi- 
bly neutron stars) in X-ray binaries (e.g. Narayan 
1996|; Es in, McClintock & Narayan p9"7| ; Menou et 
al. 1999]) to supermassive black holes in the center 



of galaxies (e.g. Narayan, Yi, & Mahadevan 1995, 
Lasota et al. |1996 , Narayan et al. 1999 ; Gammie, 
Narayan, & Blandford [jjgjg ) . 

Despite their structure being inherently two- 
dimensional, much theoretical work on ADAFs still 
relies on a vertically-integrated approach, which, al- 
though questionable, may indeed catch some of the 
essential properties of the model, as its successful 
application to various sources shows. The standard 
ADAF picture was, however, challenged in a recent 
paper by Blandford, & Begelman ( 1999 ). Since, at 
least in self-similar ADAFs, the Bernoulli number 



ried out by the wind must exceed by orders of mag- 
nitude that which is crossing the horizon, turning 
ADAFs into Advection-Dominated Inflow/Outflow 
Solutions, or ADIOS. Although the positiveness of 
the Bernoulli number is only a necessary (but not 
sufficient) condition for starting an outflow, this ar- 
gument poses a serious problem to the ADAF model 
for black hole accretion and must be addressed care- 
fully 

Goal of this Letter is to investigate how, and to 
which extent, the inclusion of the source of ADAF 
material affects the Bernoulli number and the on- 
set of a wind. Advection-dominated flows in black 
hole X-ray binaries (BHXBs) are most likely pro- 
duced by the evaporation of a Sha kura-Sunyaev disk 
(SSD, Shakura, & Sunyaev 1973 ), as observational 
and theoretical arguments suggest (e.g. Narayan, Ma- 
hadevan, & Quataert 1999|) . The evaporation pro- 
cess is not fully und erstoo d as yet (see e.g. Meyer, 
Sz Meyer-Hofmeister |1994 ; Honma 1996| ; Dullemond 
1999; Rozariska 1999| ) and definitely needs to be mod- 



is positive (as already noted by Narayan, & Yi 1994 , 
1995a), they suggested that a large outflow may form. 
If this is the case, they have shown that the mass car- 



eled in at least two spatial dimensions. Here we just 
assume a simple, analytical law for the evaporation 
rate and limit ourselves to a vertically-integrated de- 
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scription of the hot flow. A quite general argument, 
based on angular momentum conservation, indicates 
that purely inflowing solutions can not exist if the ac- 
cretion rate decreases with radius. Numerical mod- 
els, computed for several values of the a-viscosity 
parameter and of the transition radius Rq, support 
this conclusion. In all of them a stagnation radius 
(where the radial velocity vanishes) separates an in- 
ner inflowing region from an outer transsonic wind. 
We find that the Bernoulli number for the infalling 
gas is negative if the transition radius is less than 
~ 100 Schwarzschild radii. In these solutions about 
1/3 of the mass carried inwards from large radii by 
the thin disk reaches the horizon. More extended in- 
flows, with Rq ^ 100R Sl have a region of positive 
Bernoulli number, and are likely to be replaced by an 
ADIOS. 

2. THE MODEL 

We consider a highly idealized model for a hybrid 
accretion flow, in which a hot, advection-dominated 
phase coexists with a Shakura-Sunyaev disk. All the 
ADAF gas is assumed to be supplied by the evapo- 
ration of the surface layers of the SSD which extends 
down to the transition radius Rq. In the following 
we will not be concerned with the SSD anymore and 
focus our attention on the hot component. The in- 
ner rim is chosen to be at Ri n = 3GM/c 2 (M is the 
hole mass), and Ro > Ri n . The advection-dominated 
flow is described by the usual stationary, vertically- 
integrated equations (see e.g. Narayan, Kato, & 
Honma 1997] ), which now include the energy and mo- 
mentum exchange between the hot gas and the evap- 
orating material; the standard a-prescription for vis- 
cosity is retained, the pseudo-Newtonian potential is 
used to describe the gravitational field and we neglect 
radiative losses. The physics of the evaporation pro- 
cess is still unclear, so we just assume that matter is 
lost from the SSD (per unit area and time) according 
to the simple law 
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The continuity equation reads 
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(2) 



which can be immediately integrated to yield 
—2itRY,vr = Madaf(R)- In the previous expres- 
sions £ is the surface density and vr the radial veloc- 
ity, chosen to be negative for matter flowing towards 
the hole. The ADAF accretion rate now depends on 
R and is given by 



Madaf(R) 



M (R/R )-t - 



M, 



out 



R ^ Ro 

R < Rq 



(3) 

where M out = (2vr KLvr)\r^ M = 2^- 1 &qR 2 q = 
Mssd(R - » oo) is total accretion rate, and we con- 
sider only the case £ > 0, so Madaf decreases 
with R. If Mout > 0, Madaf becomes negative for 
R > R st = [1 + (M - M out )/M out ] l /tRQ and the gas 
crosses the horizon at a rate Mj n = M — M ou t < M; 
this implies that also vr has to switch sign at the 
stagnation radius R s t- 

In the following we assume that the injected mate- 
rial rotates at the Keplerian angular speed Qk an d 
that the rising gas elements move predominantly in 
the vertical direction with velocity <C QrR- Since 
the ADAF is expected to rotate at O, < Qk, the 
difference in the circular speeds produces a torque. 
Conservation of angular momentum then implies 



d(R 2 Q) d ( ,dfi\ , s „ 

(4) 

where v = 2/3ac s H is the viscosity coefficient, H 
the flow half-thickness and c s the isothermal sound 
speed. Other possible sources of friction between the 
SSD and the ADAF like, e.g., magnetic stresses, have 
been neglected. 

The ADAF has to spend part of its energy to heat 
up the injected gas which has initially the same tem- 
perature of the Shakura-Sunyaev disk. At the same 
time heat is produced by frictional dissipation. In- 
cluding both these effects, the local energy balance 
takes the form 
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where 7 is the adiabatic index and we assumed that 
the rising gas elements do not suffer any energy loss 
before thermalizing with the hot plasma. 
Finally, the radial force balance is expressed as 



2 dln(S/F) __ a 

Hr ~ 



vr; 
(6) 

the last term arises because the injected mass has 
zero radial momentum and is usually negligible. As 
expected, Eqs. @, @-@ reduce to the standard 
ADAF form for a = . 
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3. GLOBAL SOLUTIONS FOR THE HOT FLOW 

It can be easily shown that, under the usual as- 
sumptions, the flow equations admit a self-similar 
solution. However, at variance with standard 
advection-dominated models, a self-similar regime is 
not always allowed for evaporation-fed ADAFs. In 
fact, at least if < Q < f^, angular momentum con- 
servation implies (1 — 2£)/vr < 0, so inflowing self- 
similar solutions are ruled out unless £ < 1/2. The 
reason for this is as follows. If the amount of mass 
injected into the ADAF drops off too quickly with in- 
creasing R, viscous stresses can not transport enough 
angular momentum outwards because the density be- 
comes too low at large radii. The ADAF can not get 
rid of its angular momentum by transferring it to the 
SSD either, since the cold disk rotates faster. On a 
physical basis, it seems therefore unlikely that purely 
inflowing global solutions could exist. The only pos- 
sibility for the flow to transport efficiently angular 
momentum to infinity is to advect it, reversing its 
motion from a certain radius onwards. 

Eqs. (|2|), have been solved numerically us- 

ing a Henyey relaxation method (Dullemond, & Tur- 
olla 1998| ). Despite several attempts, no purely in- 
flowing solution was found for £ > 1/2, as the previ- 
ous argument predicts. In all models both the outflow 
and the inflow are transsonic. The presence of three 
critical points (the two sonic radii and the stagnation 
radius) reduces the number of the boundary condition 
from five (there are three first order and one second 
order differential equations) to two, the other three 
being replaced by regularity conditions at the critical 
points. Since the viscosity must be well-behaved at 
both the inner and outer edge (the no-torque condi- 
tion, see Narayan, Kato, & Honma |1997| ), we require 



dlog ft 
dlogR 



at R — Rin and R — R, 



out 



(7) 



This choice introduces no further degree of freedom. 
The solution depends only on the function &(R), on 
a and 7. Both the mass loss rate at large radii M ou t, 
and the stagnation radius follow from the calculation 
and are found to obey the analytical expression for 
R s t derived in §|| to high accuracy. 

We have computed several series of models with 
0.1 < a < 1, 1/2 < £ < 3/2 and 7 = 3/2, varying 
R in the range 20R S -W00R S . As Eqs. (|)-(|) show, 
vr, ft and c s are independent of &q and M; the same 
is for the ratio Mi n /M out of accretion to mass loss 
rates. The radial dependence of vr and c s for a typ- 
ical run is plotted in Fig. [I]; solutions with different 
values of a and £ show the same general behaviour. 
Only a limited £-range was considered here, but we 
have verified that steeper (e.g. exponential) injection 
laws give quite similar results. The accretion to mass 



loss ratio goes from Mi n /M out ~ 0.5 to ~ 0.3 as Rq 
increases from 20 to 500 R s and is not much depen- 
dent of a and £. Correspondingly, R st is in all cases 
about 1.5-Ro- As expected, for Rq ^ 500R S the inflow 
closely resembles standard ADAFs and, in particu- 
lar, has a nearly self-similar region at intermediate 
radii. For smaller values of Rq the accretion flow has 
no room to attain self-similarity, being squeezed in 
between the sonic and the stagnation radius. 

4. DISCUSSION 

Our solutions share with ADIOS the property that 
a significant (although very different) fraction of the 
material is expelled in a wind. The two models, how- 
ever, differ substantially in many respects. 

Blandford, & Begelman ( |1999| ) constructed self- 
similar advection-dominated solutions for which the 
Bernoulli number Be is negative assuming that 
Madaf oc R p with p > 0. They concluded that 
all the mass which does not reach the horizon es- 
capes in a wind more massive than the ADAF by 
orders of magnitude. Although no detailed model 
including the wind has been presented as yet, pre- 
liminary hydrodynamical calculations seem indeed to 
support the original suggestion that the mass inflow 
rate increases with radius (Stone, Pringle, & Begel- 
man 1999| ). In addition, 2-D simulations by Igu- 
menshchev, & Abramowicz ( |1999| ) have shown that 
ADIOS-like solutions are present for large a. In our 
model the outflow is a direct consequence of hav- 
ing assumed that the ADAF material is supplied by 
the SSD. The accreting gas has to transfer angular 
momentum to larger radii and, since Madaf is de- 
creasing with increasing R in the evaporation region, 
this can not be done by viscous stresses alone. From 
the stagnation radius onwards, the gas is centrifu- 
gally accelerated away from the hole, carrying angu- 
lar momentum with it: no stationary solution would 
be possible without an outflowing region. This is not 
related to the positiveness of the Bernoulli number. 
Moreover, the wind is not produced by the ADAF it- 
self but originates directly from the evaporating SSD 
material. 

The basic objection Blandford and Begelman raised 
to standard ADAFs is very general and concerns also 
the inflowing part of evaporation-fed models. The 
Bernoulli number for our solutions is shown in Fig. 
[||. In the wind region the Bernoulli number becomes 
positive shortly beyond the stagnation radius, while 
Be can be either positive or negative for the infalling 
gas, depending on Rq. Solutions with a small transi- 
tion radius have always Be < 0. This can be under- 
stood considering the energy the hot flow transfers 
to the cold injected material (last term in eq. ||). 
The evaporation stops at Rq but the inflowing gas 
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will stay cooler (with respect to a standard ADAF) 
a bit further down. For Rq ^ 100R S viscous dissi- 
pation has no time to heat the flow sufficiently be- 
fore it crosses the horizon, keeping Be negative (see 
also Nakamura 1995 for a discussion on the sign of 
the Bernoulli number in non-self-similar ADAF solu- 
tions). As Ro increases above ~ W0R s the Bernoulli 
number is still negative close to Ro but it flips sign as 
soon as the inflowing gas has become hot enough. An 
ADIOS-type wind may be expected where Be > 0. 

The inclusion of a source for the hot gas points 
towards the existence of a more general class of 
advection-dominated flows, with somewhat interme- 
diate characteristics between standard ADAFs and 
ADIOS. In the light of the relation to (and pos- 
sible inclusion of) ADIOS-type flows, we refer to 
these models as "Consistent Inflow And Outflows", 
or CIAOs. 

Our vertically-integrated CIAO model has to be 
supported by detailed 2-D and 3-D hydro dynamical 
calculations, but the present analysis suggests that 
for i?o <J 100i2 s the inflowing matter is gravitation- 
ally bound and no ADIOS-like wind is expected to 
be present. CIAOs can supply the central hole with 
about 1/3 to 1/4 of the mass originally carried in- 



wards by the thin disk. This may be relevant for 
applications of advection-dominated flows to BHXBs 
and AGNs. It has been often proposed (e.g. Narayan 
Esin, McClintock, 



1996 



Narayan 1997; Belloni, 



et al. 1997; Narayan, Mahadevan & Quataert 1999) 
that the different spectral states observed in BHXBs 
may be explained in terms of a bimodal disk in which 
the transition radius varies with time. The inferred 
values of Ro are in between few R s and ~ 100 R s . Re- 
cent optical/UV observations also indicate that the 
inner edge of the thin disk in the nucleus of M81 and 
NGC 4579 lies at ~ 100 R s (Quataert, et al. |1999|) . 
Although a more detailed analysis is definitely re- 
quired before any firm conclusion can be drawn, it 
is interesting to note that the maximum value of the 
transition radius for which our models have Be < is 
close to the observed limit mentioned above. CIAOs 
may comfortably provide the hot accretion flow in all 
the range of Ro implied by observations. At the same 
time, the absence of hot flows with an extent ^ 100i? s 
can be interpreted in terms of the onset of a strong 
wind which, as in ADIOS, blows off the accreting part 
of the CIAO and reduces dramatically the accretion 
rate onto the black hole. 
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Fig. 1. — The radial run of \vr\ (dotted/full line for negative/positive values) and c s (dashed line), both in units of c, for 
a — 0.5, £ = 0.75 and three values of Ro- From top to bottom Mi„/M out = 0.44, 0.37 and 0.35. The arrows mark the radial flow 
direction. 




Fig. 2. — Same as in fig. |l] for \Be\ (in units of c 2 ). Be is negative for R < R s t in the middle panel, but the dip at ~ 20R S 
signals that it is going soon to switch sign in the inflowing region. 



